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The fluctuation Conduction of NSN and SNS contacts above T c are analyzed. For NSN contacts, 
both Aslamazov - Larkin and Maki - Thompson corrections to the conduction are found to be of 
the same order and diverge for T* < T c according to the law (T — T c *) _1 . For SNS contacts, 
the Aslamazov - Larkin correction vanishes, while the Maki - Thompson correction is essential for 
contacts shorter than the phase-breaking length. 
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I. INTRODUCTION 

Recently, mesoscopic superconducting - normal metal(S/N) systems have attracted a great deal of attention 0] - 
0. In particular, it was shown that their conductance exhibits oscillatory behavior in mag netic field § - g and 
nonmonotonic temperature and voltage dependences [^) . The reason for this behaviour is the effect superconducting 
correlations have on the electrons in the normal metal. The physics of these effects is simillar to the physics of 
corrections to the conductivity resulting from superconducting fluctuations above T c |Q, |fl5|| , fl6|| . In particular, the 
authors of ref 0] have shown that the nonmonotonic temperature dependence of conductance in S /N systems is the 
result of competition between the contribution from the modified density of states and a contribution which is similar 
to the Maki - Thompson (MT) fluctuation conductivity above T c JT^j. Therefore, it is of interest to calculate the 
conductance of different S/N systems above T c taking into account superconducting fluctuations. 

Despite the large number of papers concerned with superconducting fluctuations in macroscopic samples, very few 
authors have considered superconducting fluctuations in contacts. In particular, Kulik || considered the effect of 
superconducting fluctuations on the density of states and on the current in a tunnel junction. Zaitsev ]l0| considered 
the fluctuation corrections to the conductance of very short superconducting microbridges. However, these studies 
have revealed only two types of fluctuation corrections in uniform systems; the correction due to the modified density 
of states and the MT correction, which represents the effect of fluctuational Cooper pairs on the conduction of normal 
electrons. They did not reveal the Aslamazov - Larkin (AL) correction |Q which represents the direct contribution 
of fluctuational Cooper pairs to the current. 

In this paper, we consider the effects of superconducting fluctuations on the conductance of mesoscopic NSN and 
SNS contacts of various lengths. In the case of NSN contacts, we find that the AL and MT corrections are of the 
same order of magnitude. In the case of SNS contacts, the conductance is determined by the MT correction, which 
penetrates into the contact, from the electrodes, a distance upto the phase-breaking length L v . 

II. BASIC EQUATIONS 

The expressions determining the superconducting corrections to the conductivity are easily obtained by a trivial 
extension of the Aslamazov - Larkin and Maki - Thompson equations to inhomogeneous systems. However, as many 
people are not familiar with the diagrammatic technique used by these authors, we present here a different derivation 
based on quasiclassical Green's functions of the superconductor and the self-consistency equation with a Langevin 
source [[ll| , |l0| , [Q . One can show that the results obtained with the aid of the diagrammatic technique and the 
results presented here are identical. 

The current density in a dirty superconductor is expressed by the formula 
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where g R , g A , and g F are quasiclassical matrix Green's functions of the superconductor |l5| ], Np = mpo/2Tr 2 is the 
density of states at the Fermi level, and D = lvo/3 is the diffusion coefficient of electrons. The retarded and advanced 
Green's functions g R ( A ^ describe the energy spectrum of the superconductor, while g F also contains information about 
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the electron distribution. In the case of a time-independent electrical potential, the functions g R ( A ) (e, e' , r) obey the 
equation 



D ik v &r ) + ir ieh ^ {A) ) §R(A) ~ §R{A) h e ' h ±R(A) ) = (2) 



where 

A = 



A 
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and the matrix S^'^ = ±Tt z describes the depairing. The products of the matrix quantities also imply convolutions 
over the inner frequencies. 

The order parameter A(o>, r) satisfies the self-consistency equation containing the source of condensate fluctuations 
EH: 



A(lu, r) = ^J ^Sp [{t x - T y )g F (e + u/2, e - w/2, r)] - Ajtfw, r), 



(3) 



and the correlation function of the sources of fluctuations rj is given by 



< 7?(w,r)?7*(u/,r') >= (16tvN f )- 1 6(uj + u}')6(r - r') 

<r}(u,r) V (tj',r')>=0. (4) 

This expression implies that the sources of condenstate fluctuations are S— correlated in time. The latter equality 
is the result of randomness in the phase of superconducting fluctuations above T c . 

Since the fluctuations of the order parameter are small, the retarded and advanced Green's functions of the su- 
perconductor may be represented as the sum of corresponding normal-metal Green's functions and a small additive 
proportional to A: 

g R W(e, e') = ±2nr z S(e ~ e') + f R + A \e, e'). (5) 
Substituting Eqn. into Eqn. (|J) and making use of the orthogonality condition [g RlyA ^] 2 = I [^3|, one obtains that 

f R ( A \e,e',r) = ± J dV P R ^ A \e + e', r, r') A(e - e', r'), (6) 

where the kernels P R ^ are determined by the equations 

(Dd 2 /dr 2 + ie - T) P R {e, r, r') = 2iS(r - r') 

P A (e,r,r') = -P R (-e,r,r'). (7) 

Note that the f R W matrices contain no diagonal components. The correction to the diagonal components, which 
determines the density of states, is proportional to the order parameter squared. However, this correction is small for 
the case under consideration and will be neglected by us. 

As the local electron distribution is assumed to be equilibrium, the function g F may be expressed in terms of g R 
and g A via the relationship 

g F (e 7 e<) = g R (e,e')n(e') - h(e)g A (e,e% (8) 

where 



n(e, r) 



ri(e — e<p(r)) 

n(e + e<j>(r)) J ' 



n(e) = tanh(e/2T). (9) 
and 4>(r) is the electric potential. Substituting Eqns. fl|), (||), and (|^) into the self-consistency equation (||) gives 
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A(w,r) = - / tfV / — P R (2e,r,r')n(e-oj/2 + e<j)(r)) 
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+P A (2e, r, r')n(e + w/2 - e0(r)) A(w, r') - An(u, r). 



(10) 



Consider the case where all the relevant length scales are much larger than the characteristic length £o ~ {D/Tc) 1 / 2 
and to T c . Then A(r') may be expanded in powers of r — r' to quadratic terms. Substituting the expression for 
pR{ A ) i n infinite space, 



d 3 q exp[iq(r — r')] 



> 3 e±i{Dq 2 + Y)/2 1 



(11) 



into the self-consistency equation (|l0|) , one obtains the nonstationary Ginzburg - Landau - Langevin equation in the 
form 
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D q^2 +iLJ - 2 *e0(r) - tT c - Tj A(w, r) = 16Try(w, r), 
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(12) 



where T c is the BCS transition temperature and r = (8/n)(T — T c )/T c . This equation is well known in the theory of 
superconducting fluctuations. 

Consider now the expression for the current (Q). Substituting Eqns. (|^) and (||), one obtains: 



j — jra + JAL + jMTi 



(13) 



where 
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j„ = -N F De j deSp ( f z — ) = a„E 



Oh 
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a n represents the normal-state conduction. The second term represents the regular Aslamazov - Larkin (AL) correction 
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and the third term represents the anomalous Maki - Thompson (MT) correction 
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First consider the AL correction. Substituting Eqn. (^) into Eqn. (^) and performing the averaging over the 
fluctuations of the order parameter, one obtains 



Jal = ^eN F D I ^ I <f n I d*r 2 



2-k J 2n 



P*(2e-w,r,n) 



dP R {2e-uj,r,Y 2 ) 
dr 



-P A {2e-uj,r,r 1 ) 



dP A (2e-u,r, r 2 ) 
dr 



[< A*(w, ri)A(-w, r 2 ) > n(e + e$ 



< A(a;,r 1 )A*(-w,r 2 ) > n(e 



(17) 



Equation (|T^) may be simplified, when the characteristic length scales are much larger than £oj by setting ri = r 
in the correlators in second factor of the integrand of ( |17|) and expanding in powers of r 2 — r to linear terms. Making 
use of Eqn. flU) for P R ^ in the infinite space, Eqn. (|l7|) is easily shown to be of the form 
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This is just the standard Ginzburg - Landau expression for the current. Combined with Eqn. (|12[), it gives the 
correction to the current due to fluctuations, which include the effects of a nonlinear electric field. 
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We now calculate the AL correction to linear terms in the electric field. Using the functions K r (uj, p, p'), the 
Green's function of Eqn. ( p"2| ) with zero potential, and K A (ui, p, p') = K R (—u>, p, p'), then retaining only linear terms 
in the electric field, the fluctuation of the order parameter may be written in the form 

A(w,r) = J d 3 r' K R (uj,r,r')ri(uj,r') + 2ie J d 3 v' J d 3 n K R {u>, r, n)[^(n) - <f>(r)]K R (w, n, v') V (uj, r'), 

A*(-«,r) = J d 3 r']f A ( W ,r,r'MKr')-2 ie / d 3 r' J d 3 n X A ( W , r, n)[0(n) - 0(r)]X A ( W , n, r')r?*(-^, r'). (19) 

Note that in comparison with K R (u>, p, p') defined by Aslamasov and Larkin, this quantity contains an additional 
factor mpo/l6irT c . Substituting Eqn. ( |l9| ) and the correlator of the sources of condensate fluctuations Eqn. (Q) into 
Eqn. (fL8|), one obtains the linear AL correction in the form 



j AL (r) = \e 2 DT J dto J d 3 v, J ^2 °^' ^ {Wi) - ef>(v)}K R (r u r 2 ) 

- [0(r 2 ) - 4>{v)]K A (r 2 , t{)}K a {t, r 2 ). (20) 

Note that for the AL correction, the current - field relationship is substantially nonlocal, i.e., the current density at 
a given point is determined by the electric field in its vicinity, of radius £(T) ~ ^D/{T — T c ). This is a consequence 
of the large size of fluctuational Cooper pairs near T c . 

Now we proceed to the MT correction. Substituting Eqn. (||) for /r(a) into Eqn. (16), one obtains 

j M x(r) = -ne 2 N F D^ J d 3 n J d 3 r 2 J ^ (A(w,n)A*(-a;,r 2 )> 

X / ^ pR{2e ' r ' T ^ PA ^ r < T2 ^ n{t ~ W/2 + 60(r)) - (21) 



Retaining only the terms linear in the electric field the MT correction is given by the expression 

a MT (r)=8e 2 DT J d 3 n J d 3 r 2 J d 3 r' J ^K R (u, n, v')K A (u>, r 2) r') J ^-P R (2e,v,v 1 )P A (2e,v,v 2 ) 



(22) 



Unlike the AL correction, the current-field relationship of the MT correction is purely local, as the electric field 
directly affects normal electrons rather than fluctuational Cooper pairs. 

III. NSN CONTACT 

Consider the NSN contact in the shape of a narrow channel of length L>^ and cross-sectional area S connecting 
two massive electrodes. The transverse dimensions of the channel are assumed to be much smaller than £(T). Let 
the x axis be directed along the channel. The normal-state electric potential in the contact is unperturbed by 
superconducting fluctuations and has the form = —Vx/L, where V is the voltage drop across the contact. As the 
superconducting corrections to the current essentially depend on the distance from the electrodes, the AL and MT 
corrections (Eqns. ( p0| ) , $2\)) should be calculated with the unperturbed potential and then averaged over the length 
of the contact to ensure current conservation. 

First consider the AL correction. As the transverse dimensions of the channel are small, all the relevant quantities 
may be considered to be dependent only on the longitudinal coordinate x. Introduce a system of the eigenfuctions of 
the Laplace equation 

"2 



Then the function K R entering into Eqn. (g0|) , the Green's function of Eqn. ( ^2[ ) with zero boundary conditions at 
the ends of the contact, may be represented in the form 

K*M = ~E + .^ 2 i X \ y (24) 
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where €l = D/L 2 is the Thouless energy, tl = tT c /£l, and 7 = T/cl- Performing the integration in Eqn. (E 
the coordinates and frequency and making use of the relationships 



L 



ir [m — n ) 



dx ip m (x)x(p n (x) = { (25) 

L/2, m = n, 



one obtains the AL correction in the form 



2 ran 

"7 2 2 

L> m. — m 



[l-(-l) m+n ], m^ra 
0, ra = n, 

[1 - (-l) m + n ] 2 



(26) 



m 2 -nV m 6 n (9 m + 9 n ) 



(27) 



where (9 m = 7r 2 m 2 + tl + 7. In the limit tx 3> 1 ^> 7, Eqn. ( p7j ) gives the standard AL correction for the one- 
dimensional wire 



_ ^3/2 e 2 T£ l/2 

w ^ - 2 n/2 (T-T c ) 3 /2 K ' 



(28) 



The AL correction (g7|) remains finite at T = T c owing to the finite length of the contact, the transition temperature 
decreases from the bulk value to the value 



7T 7T 

T* = T c - -r e L . 

Near T* the temperature dependence of AL correction is of the form 

e 2 TV 



(29) 



SIal = a 



a = 772 774 « 0.1026. 



(30) 



Now we proceed to the MT correction. As the quantities P R ( A ) appearing in Eqn. ( p2| ) also satisfy the zero 
boundary conditions at the ends of the contacts, they have the form 



P R( f TT V-V ^n(x)cp n (x') 

{ ' ' ' e + ie L {v 2 n 2 + 1 ) 



(31) 



Substituting Eqns. (^4|) and ( |3l| ) into Eqn. (p2|), averaging over the contact length and summing the resulting series 
one obtains 



51 mt = 2e 2 V- f i (cotW' 2 - 4 



for tl + 7 > and 



SI 



T I 7 



1 



7 1/2 7 (tl+ 7 ) 1/2 



coth(r L + 7) 



1/2 _ 



MT 



2e z Vr-r 



1 



T\ l\TL+l 



1/2 



1 : — [777 -cot|r L + 7 | 

\tl +7l 7 



1/2 



(tl+ 7 ) 1/2 
^( COt ^ 1/2 -^)} 



for tl + 7 < 0. Alternatively, Eqns.(|32j) - ( |33| ) can be written in the form 

. e 2 TV 



SImt = 4- 



y I 



(32) 



(33) 



(34) 
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For large contact lengths ti > 1, the MT correction reduces to the standard equation for the one-dimensional wire 




SI MT = y V \ - — ] . (35) 

As well as the AL correction, the MT correction remains finite at T = T c and diverges at T — T* according to the 
same law 

1 e 2 TV 

51 mt = 2^Y~T*' (36) 

In Fig. 1 we plotted the temperature dependences of the MT and AL terms in dimentionless fluctuation conductances 
Smt and Sal for different values of 7 (the ratio of the depairing rate T and the Thouless energy ej,); here Smt — 
SImt/Io , Sal = SIal/Io , lo — e 2 VT/eL- As may be seen from Fig. 1, the MT contribution dominates over the 
whole temperature range. The depairing rate T can be determined from measurements of the fluctuation conductance. 





FIG. 1. Figure 1 (a), shows the MT and the AL terms for several values of 7 (ga). Figure 1 (b), shows the ration of MT/AL 
terms for several values of 7 (ga). Figure 1 (c) is a Log-Log plot of Figure 1 (a). 
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IV. SNS CONTACT 



Consider a structure of similar geometry as the structures considered perviously but with a normal-metal channel 
and superconducting electrodes. The superconductor and normal metal are characterized by the diffusion coefficients 
D s and D n and by the phase-breaking rates T s and T n , respectively. In the case of SNS contacts, the fluctuational 
Cooper pairs generated in the superconducting electrodes can penetrate into the normal metal only a distance shorter 
than £o as the BCS coupling constant A is zero in the normal metal. However, as well as in the superconducting state 
at T < T c J! , , they can affect the conduction of normal electrons at distances much larger than the phase-breaking 
length L v . In this case, the conductance of the contact depends on the geometry of the electrodes even though their 
dimensions are much larger than the transverse dimensions of the contact. 

In view of this reasoning, the AL correction is negligible in SNS contacts. The MT correction is obtained by 
averaging Eqn. (^) over the contact length, then integrating with respect to ri, V2, and r' over the bulk of both 
electrodes, each giving an independent contribution to the conductance of the contact. Hence in Eqn. (^2|), one may 
use the expressions for K R ^ A ^ in a bulk homogeneous superconductor: 

jsR, i\ 1 \T exp[iq(r - r')] 

il ^— ' lu + i[D s q 2 + tT c + T s \ 

where O is the normalization volume. 

First consider the contribution from the left electrode. Assume that the origin coincides with the left end of 
the contact. As in the case of NSN contact, all the quantities inside the channel depend only on the longitudinal 
coordinate x. As the integral ( p2| ) is dominated by |ri| and | r"2 1 of the order of £(T), i.e., much larger than the 
transverse dimensions of the contact, the quantity P R (e, x, r), where x is the coordinate of a point inside the contact 
and r is the coordinate of a point inside the left electrode, may be represented in the form 

P R W(e,x,r) = P R{A \e,-r)^(±e,x), (38) 
where P R (e, — r) is given by Eqn. ([ll]) and function ip is the solution of the equation 

d 2 tjj 

D n-j^2 + (* e ~ T n )ip(e, x) = 
with the boundary conditions ?/>(0) = 1 and ip(L) = 0. Explicitly, it is given by the expression 

, . , sinhf«;(l — x/L)] , . . . 

<Ke, x) = L : - 1 n , k(c) = 7n - ze/e L ) 1/2 . (39) 

smh k 



With these expressions and taking into account contributions from both electrodes, Eqn. (|22|) takes the form 

o I v-^ f dio 1 f de 1 , . , 

SImt = Me 2 e L TLSV- £ / — / — —————— ^ 2e , 40 

il^J 2n lu 2 + (D s q 2 +tT c + T s ) 2 J 2tt 4e 2 + (D s q 2 + T s ) 2 



where 



*(e) = i [dx^e, xM-e, x) = ^ ~ ^ ™^ 



L J cosh(2Ki) — cos(2«2) 



k\ = Re[«(e)], k 2 = Im[«;(e)]. (41) 
Integrating with respect to frequency ui in Eqn. ( pp[ ) gives 

SImt = 32e 2 e L TLSV- V 5 — f — — — ^ *(2e). (42) 



n ^ D s q 2 + tT c + T S J 2ir 4e 2 + (D s q 2 + T s ) 2 

To be specific, consider the case where the electrodes represent a film of thickness do < £,(T), then in this case the 
sum over q may be replaced by the integral 

d 2 q 
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Introducing the dimensionless integration variable 6 = 2e/e£, one arrives at the following expression: 



Ae^TLS 
tt z D s d 



(W 



e 2 



ln 



tl + Is 
1/2 



-4- arctan I — 



2 l/2 



^sinh^ - sin 2 
cosh 6*i — cos 62 

1/2 



7r 



(43) 



Assume for simplicity that the depairing rates in the normal and superconducting metals are equal (7 S = 7„ = 7). 
First consider the case of a very short contact, 7 < Tx < 1, In this case, the integral ( (43| ) is dominated by 9 ~ tl, 
so the last factor in the integrand may be set equal to 1/3. This yields 



r t 1 e 2 

Wmt = tt In I 1 

12 L> s 



r 



sv 



(44) 



To within the factor D n /D s and a numerical coefficient, the correction to the conductivity of the contact material 
is equal to the MT conductivity of the electrodes. A similar result was previously obtained by Zaitsev for short 
ScS contacts. 

Consider now the case where the contact is of intermediate length, 7 <C 1 <C tl, one obtains with logarithmic 
accuracy 



c t 1 D n e 2 SV 

SlMT= l2D- s d^ Hlh) — - 



(45) 



This expression differs from that for the short contact in that the quantity (8/tt)(T — T c ) in the logarithm is replaced 
by the Thouless energy el- 

Lastly, consider the case of a long contact, 1 < 7 < tl. In that case, the integral (|4^ ) is dominated by 6 ^> 1, so 
the second factor in the integrand may be approximated by (28) ~ 1 ^ 2 , yielding 



SV 



XT 1Dn e -1/2 

Shir = 7;TT1~ 7 r 
2 Us do T Li 



(46) 



The physical meaning of this result is that the superconducting correlations that result in the MT correction 
penetrate into the contact over the length L v ~ 7 _1 / 2 L <c L. 

Note that unlike the case of NSN contact, the correction to the total current is proportional to the cross-sectional 
area of the contact. 



V. CONCLUSIONS 



The fluctuation conductivities of NSN and SNS contacts of arbitrary lengths have been calculated. We have 
established that the fluctuation conductivity in NSN contacts consists of contributions from both the Maki - Thompson 
and Aslamasov - Larkin terms. The MT contribution dominates over the whole temperature range. However, near 
the renormalized critical temperature T* the ratio of the MT and AL terms does not contain any parameters, and 
equals about 1.56 . 

In SNS contacts the AL contribution is absent. The increase in the conductivity due to fluctuations is caused by 
the anomalous MT term containing the product of retarded and advanced functions (see Eqn. (p2|)). 

The superconducting fluctuations modiify the density of states and decrease the conductivity. In the case of NSN 
and SNS contacts analyzed by us, this decrease is small (i.e. it does not diverge as T tends to T*). However the 
correction to the conductivity due to the decrease of the density of states is essential in SNS contacts at T < T* , this 
leads to the reentrant behaviour of the conductance Q, Q, ||, and is also essential in tunnel SIS junctions and in 
layered superconductors O], |17 
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